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A  central  question  in  the  field  of  graphene-related  research  is 
how  graphene  behaves  when  it  is  patterned  at  the  nanometre 
scale  with  different  edge  geometries.  A  fundamental  shape 
relevant  to  this  question  is  the  graphene  nanoribbon  (GNR), 
a  narrow  strip  of  graphene  that  can  have  different  chirality 
depending  on  the  angle  at  which  it  is  cut.  Such  GNRs  have 
been  predicted  to  exhibit  a  wide  range  of  behaviour,  including 
tunable  energy  gaps1,2  and  the  presence  of  one-dimensional 
(ID)  edge  states3-5  with  unusual  magnetic  structure6,7.  Most 
GNRs  measured  up  to  now  have  been  characterized  by  means  of 
their  electrical  conductivity,  leaving  the  relationship  between 
electronic  structure  and  local  atomic  geometry  unclear8-10. 
Here  we  present  a  sub-nanometre-resolved  scanning  tun¬ 
nelling  microscopy  (STM)  and  spectroscopy  (STS)  study  of 
GNRs  that  allows  us  to  examine  how  GNR  electronic  struc¬ 
ture  depends  on  the  chirality  of  atomically  well-defined  GNR 
edges.  The  GNRs  used  here  were  chemically  synthesized  using 
carbon  nanotube  (CNT)  unzipping  methods  that  allow  flexible 
variation  of  GNR  width,  length,  chirality,  and  substrate11,12.  Our 
STS  measurements  reveal  the  presence  of  ID  GNR  edge  states, 
the  behaviour  of  which  matches  theoretical  expectations  for 
GNRs  of  similar  width  and  chirality,  including  width-dependent 
energy  splitting  of  the  GNR  edge  state. 

The  chirality  of  a  GNR  is  characterized  by  a  chiral  vector 
( n,m )  or,  equivalently,  by  chiral  angle  9,  as  shown  in  Fig.  la. 
GNRs  having  different  widths  and  chiralities  were  deposited  on 
a  clean  Au(lll)  surface  and  measured  using  STM.  Figure  lb 
shows  a  room  temperature  image  of  a  single  monolayer  GNR 
(GNR  height  is  determined  from  linescans,  such  as  that  shown 
in  Fig.  lb  inset;  some  multilayer  GNRs  were  observed,  but  we 
focus  here  on  monolayer  GNRs).  The  GNR  of  Fig.  lb  has  a  width 
of  23.1  nm,  a  length  greater  than  600  nm,  and  exhibits  straight, 
atomically  smooth  edges  (the  highest  quality  GNR  edges,  such  as 
those  shown  in  Figs  1  and  2,  were  observed  in  GNRs  synthesized 
as  in  ref.  11).  Such  GNRs  are  seen  to  have  a  ‘bright  stripe’ 
running  along  each  edge. 

This  stripe  marks  a  region  of  curvature  near  the  terminal 
edge  of  the  GNR  that  has  a  maximum  extension  of  ~3  A  above 
the  mid-plane  terrace  of  the  GNR  and  a  width  of  ~30A  (see 
line  scan  in  Fig.  lb  inset).  Such  edge-curvature  was  observed 
for  all  high-quality  GNRs  examined  in  this  study  (more  than 
150,  including  GNRs  deposited  onto  a  Ru(0001)  surface).  This 


Figure  1  |  Topography  of  graphene  nanoribbons  (GNRs)  on  Au(111).  a,  A 

schematic  drawing  of  an  (8, 1)  GNR.  The  chiral  vector  ( n,m )  connecting 
crystallographically  equivalent  sites  along  the  edge  defines  the  edge 
orientation  of  the  GNR  (black  arrow).  The  blue  and  red  arrows  are  the 
projections  of  the  (8, 1)  vector  onto  the  basis  vectors  of  the  graphene 
lattice.  Zigzag  and  armchair  edges  have  corresponding  chiral  angles  of 
6  =  0°  and  9  =  30°,  respectively,  whereas  the  (8, 1)  edge  has  an  chiral 
angle  of  9  =  5.8°.  b,  Constant-current  STM  image  of  a  monolayer  GNR  on 
Au(111)  at  room  temperature  (Vs  =1.5  V,  /  =  100  pA).  Inset  shows  the 
indicated  line  profile,  c,  Higher  resolution  STM  image  of  a  GNR  at  T=  7  K 
(Vs  =  0.2  V,  /  =  30  pA,  greyscale  height  map). 

is  reminiscent  of  curved  edge  structures  observed  previously  near 
graphite  step-edges13.  We  rule  out  that  these  GNRs  are  collapsed 
nanotubes  by  virtue  of  the  measured  ratio  (observed  to  be  Jt)  of 
GNR  width  to  nanotube  height  for  partially  unzipped  CNTs.  We 
further  rule  out  that  the  curved  GNR  edges  observed  here  are 
folded  graphene  boundaries  by  means  of  a  detailed  comparison  of 
terminal  curved  edges  and  actual  folded  edges  (see  Supplementary 
Information).  Low-temperature  STM  images  (Figs  lc  and  2a)  show 
finer  structure  in  both  the  interior  GNR  terrace  and  the  edge 
region.  Figure  2a,  for  example,  shows  the  atomically-resolved  edge 
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Figure  2  |  Edge  states  of  GNRs.  a,  Atomically-resolved  topography  of  the  terminal  edge  of  an  (8, 1)  GNR  with  measured  width  of  19.5  ±0.4  nm 
(Vs  =  0.3  V,  /  =  60  pA,  T  =  7  K).  b,  Structural  model  of  the  (8, 1)  GNR  edge  shown  in  a.  c,  dl/dV  spectra  of  the  GNR  edge  shown  in  a,  measured  at  different 
points  (black  dots,  as  shown)  along  a  line  perpendicular  to  the  GNR  edge  at  T=  7  K.  Inset  shows  a  higher  resolution  d//d  V  spectrum  for  the  edge  of  a  (5,  2) 
GNR  with  width  of  15.6 ±0.1  nm  (initial  tunnelling  parameters  Vs  =  0.15  V,  /=  50  pA;  wiggle  voltage  Vr.m.s.  =  2  mV).  The  dashed  lines  are  guides  to  the 
eye.  d,  dl/dV  spectra  measured  at  different  points  (red  dots,  as  shown)  along  a  line  parallel  to  the  GNR  edge  shown  in  a  at  T=  7  K  (initial  tunnelling 
parameters  for  c  and  d  are  Vs  =  0.3  V,  /  =  50  pA;  wiggle  voltage  Vr.m.s.  =  5  mV). 


region  of  a  monolayer  GNR  and  clearly  exhibits  how  the  periodic 
graphene  sheet  of  the  GNR  terminates  cleanly  and  with  atomic 
order  at  the  gold  surface. 

Such  high-resolution  images  allow  us  to  experimentally  deter¬ 
mine  the  chirality  of  GNRs,  and  to  create  structural  models  of 
observed  edge  regions.  In  Fig.  2a,  for  example,  we  see  rows  of 
protrusions  (with  a  spacing,  ~2.5A,  equivalent  to  the  graphene 
lattice)  near  the  edge  of  a  GNR  having  a  width  of  19.5  ±  0.4  nm. 
By  comparing  this  row  orientation  with  the  orientation  of  the  GNR 
terminal  edge  we  are  able  to  extract  the  GNR  chirality  (details  in 
Supplementary  Information).  The  GNR  shown  in  Fig.  2a  has  an 
(8,  1)  chirality  (equivalent  to  9  =  5.8°),  and  the  resulting  structural 
model  for  this  GNR  is  shown  in  Fig.  2b.  We  find  the  distribution  of 
GNR  chiralities  to  be  essentially  random.  This  is  consistent  with  our 
structural  data,  which  indicates  that  the  CNT  unzipping  direction 
is  very  close  to  the  axial  direction  of  the  precursor  CNTs  (see 


Supplementary  Information),  as  well  as  the  fact  that  the  precursor 
CNTs  have  a  broad  chirality  distribution14. 

We  explored  the  local  electronic  structure  of  GNR  edges  using 
STS,  in  which  AI/AV  measurement  reflects  the  energy-resolved 
local  density  of  states  (LDOS)  of  a  GNR.  Figure  2c,d  shows  d I /dV 
spectra  obtained  at  different  positions  (as  marked)  near  the  edge 
of  the  (8,  1)  GNR  pictured  in  Fig.  2a.  dl/dV  spectra  measured 
within  24  A  of  the  GNR  edge  typically  show  a  broad  gap-like  feature 
having  an  energy  width  of  ~130meV.  This  is  very  similar  to  the 
behaviour  observed  in  the  middle  of  large-scale  graphene  sheets, 
and  is  attributed  to  the  onset  of  phonon-assisted  inelastic  electron 
tunnelling15  for  |£|  >  65  meV.  This  feature  disappears  further  into 
the  interior  of  the  GNR,  as  expected,  because  of  increased  tunnelling 
to  the  Au  substrate16.  Very  close  to  the  GNR  edge,  however,  we 
observe  further  features  in  the  spectra.  The  most  dominant  of  these 
features  are  two  peaks  that  rise  up  within  the  elastic  tunnelling 
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region  (that  is  at  energies  below  the  phonon-assisted  inelastic  onset) 
and  which  straddle  zero  bias.  For  the  GNR  shown  in  Fig.  2a  (which 
has  a  width  of  19.5  ±  0.4  nm)  the  two  peaks  are  separated  in  energy 
by  a  splitting  of  A  =  23.8  ±3.2  meV.  Similar  energy-split  edge-state 
peaks  have  been  observed  in  all  the  clean  chiral  GNRs  that  we 
investigated  spectroscopically  at  low  temperature.  For  example, 
the  inset  to  Fig.  2c  shows  a  higher  resolution  spectrum  exhibiting 
energy-split  edge-state  peaks  for  a  (5,  2)  GNR  having  a  width  of 
15.6  nm  and  an  energy  splitting  of  A  =  27.6  ±  1.0  meV.  The  two 
edge-state  peaks  are  often  asymmetric  in  intensity  (depending  on 
specific  location  in  the  GNR  edge  region),  and  their  mid-point  is 
often  slightly  offset  from  Vs  =  0  (within  a  range  of  ±20meV).  As 
seen  in  the  spectra  of  Fig.  2c,  the  amplitude  of  the  peaks  grows 
as  one  moves  closer  to  the  terminal  edge  of  the  GNR,  before 
falling  abruptly  to  zero  as  the  carbon/gold  terminus  is  crossed.  The 
spatial  dependence  of  the  edge-state  peak  amplitude  as  one  moves 
perpendicular  from  the  GNR  edge  is  plotted  in  Fig.  3a  and  shows 
exponential  behaviour.  The  edge-state  spectra  also  vary  as  one 
moves  parallel  to  the  GNR  edge,  as  shown  in  Fig.  2d.  The  parallel 
dependence  of  the  edge-state  peak  amplitude  is  plotted  in  Fig.  3b, 
and  oscillates  with  an  approximate  20  A  period,  corresponding 
closely  to  the  2 1 A  periodicity  of  an  ( 8,  1 )  edge. 

We  have  also  characterized  monolayer  GNRs  having  different 
chiralities  and  widths  (the  lengths  of  the  GNRs  used  in  these 
measurements  are  greater  than  500  nm).  In  Fig.  3c,  we  plot  the 
width  dependence  of  the  measured  energy  gap  of  GNR  edge  states 
for  a  broad  range  of  chirality  (3.7°  <  0  <  16.1°).  The  measured 
edge-state  energy  splitting  shows  a  clear  inverse  correlation  with 
GNR  width.  Our  gap  values  tend  to  be  smaller  than  those 
observed  previously  for  lithographically  patterned  GNRs  (probably 
because  of  uncertainty  in  the  edge  structure  of  lithographically 
obtained  GNRs;  ref.  9). 

The  high  quality  of  the  atomically  well-defined  edge  structures 
observed  here  allows  us  to  quantitatively  compare  our  experimental 
data  to  theoretical  calculations  of  the  electronic  structure  of 
chiral  GNRs.  We  find  that  the  spectroscopic  features  we  observe 
correspond  closely  to  the  spatial-  and  energy-dependence  predicted 
for  ID  spin-polarized  edge  states  coupled  across  the  width  of  a 
chiral  GNR.  This  behaviour  is  quite  different  from  the  properties 
observed  previously  for  graphite  step  edges,  armchair  nanoribbons, 
and  comparatively  less  ordered  graphene  platelet  edges,  where  no 
magnetism-induced  energy  splitting  has  been  seen17-20. 

To  compare  our  experimental  data  with  theoretical  predictions 
for  GNRs,  we  used  a  Flubbard  model  Hamiltonian,  solved  self- 
consistently  in  the  mean-field  approximation6,  for  an  (8,  1)  GNR 
having  the  same  width  as  the  actual  (8,  1)  GNR  shown  in  Fig.  2a. 
The  Hamiltonian: 

H  =  -ty^[CjJtCj»  +  h.c.]+I/y^  n;t  «4  ( 1 ) 

(y).tf 

consists  of  a  one-orbital  nearest-neighbour  tight-binding  Hamilto¬ 
nian  with  an  on-site  Coulomb  repulsion  term.  In  this  expression  ci(Tf 
and  Cja  are  operators  that  create  and  annihilate  an  electron  with  spin 
a  at  the  nearest-neighbour  sites  i  and  j  respectively,  t  =  2.7  eV  is  a 
hopping  integral21,  nia  =  c,„ 1  c,„  is  the  spin-resolved  electron  density 
at  site  i,  and  U  is  an  on-site  Coulomb  repulsion.  This  GNR  model 
is  defined  only  by  the  Jt-bonding  network.  The  terminal  a  -bonds  at 
the  GNR  edges  are  considered  to  be  passivated  and  do  not  alter  the 
JT-system  (this  should,  in  general,  correctly  model  a  range  of  differ¬ 
ent  possible  edge-adsorbate  bonding  configurations7’22,  including 
the  likely  oxygen-related  functional  group  termination  of  our  own 
GNRs;  ref.  11).  The  out-of- plane  curvature  seen  experimentally 
near  GNR  edges  is  not  included  in  this  model  because  the  measured 
radii  of  curvature  are  sufficiently  large  (>20  A)  that  they  are  not 
expected  to  significantly  affect  GNR  electronic  structure23  (we 
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Figure  3  |  Position-  and  width-dependent  edge-state  properties,  a,  Solid 
blue  dots  show  experimental  edge-state  peak  amplitude  at  points  along  a 
line  perpendicular  to  the  carbon/gold  edge  terminus  (same  positions  as 
shown  in  Fig.  2c).  Peak  amplitude  and  energies  were  determined  by  fitting 
Lorentzian  curves  to  the  two  peaks  observed  in  the  measured  spectra  at 
each  location  in  Fig.  2c  over  the  range  —30  mV  <  Vs  <  30  mV.  The  energy 
positions  of  these  peaks  were  found  to  be  6.7  ±1.6  mV  and 
—17.2  ±2.2  mV.  The  positional  dependence  of  the  average  peak  amplitude 
of  these  two  peaks  is  plotted.  Error  bars  (shown  when  larger  than  plotted 
points)  reflect  the  range  of  Lorentzian  parameters  that  result  in  a  good  fit  to 
the  data.  Dashed  red  line  shows  the  calculated  LDOS  at  locations  spaced 
perpendicular  to  the  edge  terminus  for  an  (8, 1)  GNR  (see  text)  at  the 
energy  of  the  DOS  peak  nearest  the  band-edge.  Theoretical  LDOS  values 
include  a  single  globally  constant  offset  to  model  the  added  contribution 
from  Au  surface  LDOS,  and  a  single  globally  constant  multiplicative  factor 
to  model  the  unknown  total  area  of  the  STM  tunnel  junction,  b,  Solid  blue 
dots  show  experimental  average  edge-state  peak  amplitude  (determined 
as  in  a)  at  locations  spaced  along  a  line  parallel  to  the  carbon/gold  edge 
terminus  (same  positions  as  shown  in  Fig.  2d).  Dashed  red  line  shows  the 
theoretical  edge-state  LDOS  for  an  (8, 1)  GNR  at  points  parallel  to  the  edge 
terminus  (calculated  as  in  a).  The  edge-state  LDOS  amplitude  oscillates 
parallel  to  the  edge  with  a  21  A  period,  c,  Width  dependence  of  the 
edge-state  energy  gap  of  chiral  GNRs.  From  left  to  right,  the  chiralities  of 
experimentally  measured  GNRs  are  (13, 1),  (3, 1),  (4, 1),  (5,  2),  and  (8, 1) 
respectively,  corresponding  to  a  range  of  chiral  angle  3.7°  <9  <  16.1°. 
Energy  gaps  determined  by  Lorentzian  fits  to  d//dV  peaks  (centre-to-centre 
width),  error  bars  reflect  standard  deviation  due  to  spatial  variation  in 
spectra.  GNR  width  measured  as  distance  between  the  GNR  edge 
mid-heights  on  opposite  sides,  error  bars  reflect  standard  deviation  due  to 
spatial  variation  along  GNR  axis.  The  pink  shaded  area  shows  the  predicted 
range  of  edge-state  bandgaps  as  a  function  of  width,  evaluated  for  chiral 
angles  in  the  range  0°  <  9  <  15°  (U=  0.5 1,  f  =  2.7  eV;  ref.  27). 
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Figure  4  |  Theoretical  band  structure  and  density  of  states  (DOS)  of  a  20-nm-wide  (8, 1)  GNR.  a,  Dashed  blue  line  shows  the  calculated  GNR  electronic 
structure  in  the  absence  of  electron-electron  interactions  (L/=  0).  Solid  red  line  shows  the  calculated  GNR  electronic  structure  for  U=  0.5 1  (t  =  2.7  eV). 
Finite  U  >  0  splits  degenerate  edge  states  at  £=  0  into  spin-polarized  bands,  opening  a  bandgap  (arrows),  b,  Dashed  blue  line  shows  the  (8, 1)  GNR  DOS 
for  the  U  =  0  case  in  a.  The  peak  at  £=  0  is  due  to  the  degeneracy  of  edge  states  in  the  absence  of  electron-electron  interactions.  Solid  red  line  shows  the 
(8, 1)  GNR  DOS  for  U=  0.5 1.  The  opening  of  the  bandgap  (arrows)  reflects  the  predicted  energy  splitting  due  to  the  onset  of  magnetism  in  spin-polarized 
edge  states  for  U  >  0,  and  compares  favourably  with  the  experimental  data  for  the  (8, 1)  GNR  of  Fig.  2. 


tested  this  conjecture  by  including  the  observed  curvature  in  some 
calculations,  and  found  that  it  has  no  significant  effect — either 
from  cr-Jt  coupling  or  from  pseudofield  effects — on  the  calculated 
GNR  electronic  structure).  The  effect  of  the  gold  substrate  here  is 
taken  only  as  a  charge  reservoir  that  can  slightly  shift  the  location 
of  Er  within  the  GNR  band  structure  and  reduce  the  magnitude 
of  the  effective  U  parameter  by  means  of  electrostatic  screening 
(the  experimental  charge-induced  energy  shifts  seen  here  are  within 
the  range  of  charge-induced  energy  shifts  observed  previously  for 
CNTs  on  Au;  ref.  24). 

We  first  calculated  the  GNR  electronic  structure  for  (7  =  0,  which 
effectively  omits  the  electron-electron  interactions  responsible  for 
the  onset  of  magnetic  correlations.  This  results  in  the  theoretical 
band  structure  and  density  of  states  (DOS)  shown  in  Fig.  4a, 
b  (blue  dashed  lines).  The  finite  width  of  the  GNR  leads  to  a 
family  of  subbands  in  the  band  structure,  with  no  actual  bandgap 
(Fig.  4a).  A  flat  band  at  E  =  0  due  to  localized  edge  states  leads 
to  a  strong  van  Hove  singularity  (that  is,  a  peak)  in  the  DOS  at 
E  =  0  (Fig.  4b).  The  DOS  in  this  case  does  not  resemble  what  is 
seen  experimentally.  We  next  calculated  the  (8,  1)  GNR  electronic 
structure  for  U  >  0.  Here  the  electron-electron  interactions  lift  the 
degeneracy  of  the  edge  states  by  causing  ferromagnetic  correlations 
to  develop  along  the  GNR  edges  and  antiferromagnetic  correlations 
to  develop  across  the  GNR.  This  leads  to  a  spin-polarization  of 
the  edge  states  that  splits  the  single  low-energy  peak  seen  in 
the  [7  =  0  DOS  into  a  series  of  van  Hove  singularities,  thus 
opening  up  a  gap  at  E  =  0.  Such  behaviour  is  seen  in  the  band 
structure  and  DOS  of  Fig.  4a,  b  (solid  red  lines).  We  identify 
the  lowest-energy  pair  of  van  Hove  singularities  with  the  pair  of 
peaks  observed  experimentally  near  zero  bias  for  GNR  edges.  We 
focus  our  experiment/theory  comparison  on  the  low-energy  regime 
(|£|  <  65meV)  because  higher  energy  experimental  features  are 
complicated  by  the  onset  of  phonon -assisted  inelastic  tunnelling25 
(the  low-energy  edge-state  peaks,  by  contrast,  do  not  have  the 
characteristics  of  inelastic  modes). 

We  find  that  our  experimental  spectroscopic  edge-state  data 
for  the  (8,  1)  GNR  is  in  agreement  with  model  Hamiltonian 
calculations  for  U  =  0.5 1.  The  theoretical  bandgap  of  29  meV  is  very 
close  to  the  experimentally  observed  value  of  23.8  ±  3.2  meV  (the 
value  of  [7  used  here  is  lower  than  a  value  obtained  previously  from 
a  first-principles  calculation26,  presumably  because  of  screening 


from  the  gold  substrate).  Our  experimentally  observed  energy-split 
spectroscopic  peaks  thus  provide  evidence  for  the  formation  of 
spin-polarized  edge  states  in  pristine  GNRs  (such  splitting  does 
not  arise  for  the  non-magnetic  (7  =  0  case  described  above). 
We  are  further  able  to  compare  the  spatial  dependence  of  the 
calculated  edge  states  with  the  experimentally  measured  STS  results. 
The  dashed  line  in  Fig.  3a  shows  the  theoretical  LDOS  calculated 
at  the  energy  of  the  low-energy  edge-state  peaks  as  one  moves 
perpendicularly  away  from  the  GNR  edge  and  into  the  (8,  1)  GNR 
interior.  The  predicted  exponential  decay  length  of  ~  12  A  is  in 
reasonable  agreement  with  the  experimental  data.  The  variation 
seen  in  the  calculated  LDOS  of  the  edge  state  in  the  direction  parallel 
to  the  GNR  edge  also  compares  favourably  with  our  experimental 
observations  (Fig.  3b).  The  oscillation  in  edge-state  amplitude  is 
seen  to  arise  from  ‘kinks’  in  the  zigzag  edge  structure  resulting  from 
the  chiral  nature  of  the  (8,  1)  GNR  edge  (see  the  edge  structure  of 
Fig.  2b,  each  of  the  two  dips  in  spectroscopic  amplitude  occurs  at 
the  location  of  a  kink). 

We  are  similarly  able  to  compare  the  GNR  width  dependence 
of  our  experimentally  measured  edge-state  gaps  to  theoretical 
calculations.  As  the  measured  GNRs  having  different  widths 
also  have  different  chiralities,  we  have  calculated  the  theoretical 
edge-state  gap  versus  width  behaviour  over  the  chirality  range 
0°  <  9  <  15°  (it  will  be  useful  in  the  future  to  measure  the 
local  electronic  properties  of  armchair  GNRs  (9  =  30°),  which  are 
predicted  to  have  no  edge  states).  The  pink  shaded  region  in  Fig.  3c 
shows  the  results  of  our  calculations,  and  compares  favourably 
with  our  experimentally  observed  width-dependent  edge-state  gap. 
This  provides  strong  evidence  that  the  edge-state  gap  we  observe 
experimentally  is  not  a  local  effect,  as  might  occur,  say,  in  response 
to  some  unknown  molecules  bound  to  the  GNR  edge,  but  rather 
depends  on  the  full  GNR  electronic  structure,  including  interaction 
between  the  edges. 

Methods 

The  GNRs  shown  in  this  study  were  produced  by  unzipping  carbon  nanotubes11 . 
GNRs  were  deposited  onto  clean  Au(lll)  surfaces  using  a  spin-coating  method. 
Au(lll)  substrates  were  first  cleaned  by  sputtering  and  annealing  in  ultra-high 
vacuum  (UHV)  before  spin-coating.  The  samples  were  then  transferred  into  the 
UHV  chamber  of  our  STM  system  (base  pressure  ~2.0  x  10~10  torr).  After  heat 
treatment  up  to  500  °C  in  UHV,  the  samples  were  directly  transferred  onto  the 
STM  stage  in  the  same  chamber  for  measurements. 
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STM  measurements  were  performed  using  a  home-built  STM  held  at  low 
temperature  (T  =  7  K)  for  maximum  spatial  and  spectroscopic  resolution.  STM 
topography  was  obtained  in  constant-current  mode  using  a  Ptlr  tip,  and  dl/dV 
spectra  were  measured  using  lock-in  detection  of  the  a.c.  tunnelling  current  driven 
by  a  451  Hz,  1-5  mV  (r.m.s.)  signal  added  to  the  junction  bias  (the  sample  potential 
referenced  to  the  tip)  under  open-loop  conditions.  We  also  performed  large-scale 
topographic  surveys  of  GNR  samples  before  low-temperature  measurement.  This 
was  done  in  an  Omicron  variable  temperature  STM  in  UHV  at  room  temperature. 
The  three  large-scale  topographic  STM  images  shown  (Fig.  lb  and  Supplementary 
Figs  Sla,b)  were  obtained  in  the  Omicron  STM  at  room  temperature. 
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